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After delineating the physical regimes which vortex lattices encounter in rotating Bose-Einstein 
condensates as the rotation rate, Q., increases, we derive the normal modes of the vortex lattice 
in two dimensions at zero temperature. Taking into account effects of the finite compressibility, 
we find an inertial mode of frequency > 20., and a primarily transverse Tkachenko mode, whose 
frequency goes from being linear in the wave vector in the slowly rotating regime, where is small 
compared with the lowest compressional mode frequency, to quadratic in the wave vector in the 
opposite limit. We calculate the correlation functions of vortex displacements and phase, density 
and superfluid velocities, and find that the zero-point excitations of the soft quadratic Tkachenko 
modes lead in a large system to a loss of long range phase correlations, growing logarithmically 
with distance, and hence lead to a fragmented state at zero temperature. The vortex positional 
ordering is preserved at zero temperature, but the thermally excited Tkachenko modes cause the 
relative positional fluctuations to grow logarithmically with separation at finite temperature. The 
superfluid density, defined in terms of the transverse velocity autocorrelation function, vanishes at 
all temperatures. Finally we construct the long wavelength single particle Green's function in the 
rotating system and calculate the condensate depletion as a function of temperature. 



I. INTRODUCTION 

Under rapid rotation, a superfluid forms a triangular lattice of quantized vortices carrying the angular momentum of 
the system. Such structure was seen experimentally in superfluid helium in Refs. [1]. Rapidly rotating Bose-Einstein 
condensates of cold atoms [2-4] open up the study of the physics of vortex lattices in regimes well beyond those 
achievable in superfluid helium [5-7]. With increasing rotational frequency, fi, condensates go through a number of 
different regimes. At VL small compared with the lowest compressional frequencies, U, <ti sko, where s is the sound 
velocity and fco ~ t^/R± is the lowest wavenumber in the finite geometry, with Rj_ the size of the system perpendicular 
to the rotation axis, the system is in the "stiff' Thomas- Fermi regime, and responds to rotation effectively as an 
incompressible fluid. As fl becomes larger than the lowest sound mode frequencies, fl ^ sko, or essentially that 
the outer edge of the cloud moves supersonically, the system enters the "soft" Thomas- Fermi regime, where the 
compressibility becomes important in the response to rotation. 

When, in harmonically trapped condensates, fl approaches the transverse trapping frequency, uj±, the centrifugal 
force begins to balance the trapping force, and the system flattens towards a lower density and therefore more 
compressible, two-dimensional cloud. For f2 ^ ms, where m is the particle mass, the condensate wave function 
becomes formed primarily of particle orbits in the lowest Landau level of the Coriolis force - the mean flcld quantum 
Hall state [8] - with approximate order parameter, for rotation about the z axis, 

*(r) =Cnj(M- Uj)e"l"l'/^''i, u^x + iy, (1) 

where the uj — xj + iyj are the vortex locations in the usual complex notation, d±^ = {mujj_)^^^^ is the transverse 
oscillator length, and C is a normalization constant. The structure of the trapped condensate in the axial direction 
approaches a Gaussian shape as \Omega uj±. Although Eq. (1) predicts that the transverse structure is also 
Gaussian, the transverse structure, if Thomas-Fermi in the non-rotating cloud, remains Thomas-Fermi in this limit 
[9], owing to a small admixture of higher Landau levels in the order parameter. Experiments in this regimes are 
reported in Ref. [6]. With further increase of O, the vortex lattice is expected to melt, at the point where the number 
of vortices becomes of order ten percent of the number of particles, [10,11], and the system eventually enters new 
highly correlated bosonic quantum Hall many-particle states, no longer describable in mean held [12-16]. 

Prior to discussing the modes of the lattice, it is useful to lay out the demarcations between the various regimes in 
a harmonically trapped gas. We assume that the interaction between the particles, of total number N, is described 
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by a repulsive s-wave interaction parameter g = 47ras/TO, where as is the scattering length; we use units in which 
h = 1. In Thomas- Fermi, the transverse radius, i?_L, is given by [17], 



where t = [(15A^&as/(ij^)((jj2/ti)^)]^/^, and lj^ is the axial trapping frequency. The factor h, > 1, describes the 
renormalization of the interaction energy of long wavelength density fluctuations in the system [18]. Furthermore, 

" ' (3) 



ms(0)^ T u)±_ 

where s(0) is the sound velocity in the center of the cloud. If we write fco — a/ Rl, where a ~ 5.45 [19], the criterion 
for being in the soft regime, that fl/sko be large, becomes 

""^"^ » ^ - 0.13. (4) 

The experiments of Ref [5] on Tkachenko modes reach Cl/sko ^ 1.15. 

The criterion to be in the mean field quantum Hall regime, Q./ms^ ^ 1, can be conveniently written in terms of 
the filling factor, v = N/Ny, where Ny is the total number of vortices in the rotating cloud. Since 



where n„ = m^/n is the density of vortices, we find 



(0)2 j/2/3 \uj^) \lbbasuj. 



(6) 



This result is independent of the number of particles, as long as the system is in the Thomas-Fermi regime. For the 
parameters u}^_/u!z = 8.3/5.2 of the Tkachenko mode experiment of Ref. [5] in ^'^Rb, we find f2/ms(0)2 ~ 200/i^2/3 jj^ 
the limit Vt u}±. 

The vortex lattice supports a number of modes, first discussed by Tkachenko [20] for a two-dimensional incom- 
pressible fiuid, and later in Ref. [21] at finite temperature with full effects of the normal fluid, dissipation and Kelvin 
oscillations of the vortex lines in three dimensions. The low frequency in-plane Tkachenko mode is an elliptically po- 
larized oscillation of the vortex lines, with the semi-major axis of the ellipse orthogonal to the direction of propagation. 
The Tkachenko mode is linear at small wave vector, k, in the transverse plane. 
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where n is the particle density, and C2 is the elastic shear modulus of the vortex lattice; at slow rotation, C2 = nfl/8. 
In the soft regime, the dispersion relation instead becomes quadratic [11,22], 

Were it possible to rotate helium sufficiently rapidly one would also see this very long wavelength quadratic behavior. 

Such Tkachenko soft modes can play havoc with the stability of a large system, causing loss of long range phase 
coherence even at zero temperature; they are eventually responsible for the melting of the lattice [11]. In a recent 
paper [22] we derived the modes of the vortex lattice for all rotation rates, through constructing the conservation 
laws and superfluid acceleration equation describing the long wavelength behavior of the system. In this paper we 
focus on deriving the correlation functions of density, superfiuid phase and velocity, and vortex displacements from 
equilibrium, which enable us to understand the effects of the soft infrared structure on the stability of the superfiuid 
and lattice. This work is a generalization of Ref. [23], which discussed the effects of the oscillations of the vortex 
lines at finite temperature in liquid helium on the long ranged phase correlations of the superfluid. As we shall see, 
the long wavelength Tkachenko modes lead to fragmentation of the condensate, even at zero temperature. Whether 
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the system loses phase coherence over its volume or the lattice melts first depends on the number of particles in the 
system and its rotation rate. 

In Sec. II, we review the basic equations describing the dynamics, restricting the analysis to linearized motion in 
two dimensions, and neglecting the normal component of the supcrflnid as well as dissipative terms. The analysis of 
the full three-dimensional problem will be published separately [24]. In Sec. Ill we construct the correlation functions 
of the physical quantities of interest, and in Sec. IV study the condensate depletion by constructing the single particle 
Green's function for the rotating superfluid. 



II. CONDENSATE PHASE AND CONSERVATION LAWS 



Let us first recapitulate the conservation laws and equation for the superfluid phase which govern the long wavelength 

behavior of the system. While much of this material has been given in Rcf. [22] , we include it here in order to facilitate 
the derivation of the correlation functions. The basic formalism given here applies to general bosonic superfluids. 



A. Condensate phsise 

We work in the frame co-rotating with the lattice, and describe the deviations of the vortices from their home 
positions by the continuum displacement field, e{r,t). In linear order in the vortex displacements, the long wave- 
length superfluid velocity, v{r,t), can be written, following Ref. [23], in terms of the long wavelength vortex-lattice 
displacement field and the phase $(r, t) of the order parameter, as 

v + 2nx e = V^/m; (9) 

The curl of this equation is 

Vxv = -2n'V-e. (10) 
The origin of Eq. (9) is the law of conservation of vorticity, w = V x v: 

^+Vx(^xe) = 0, (11) 

where here e tells the rate at which the vorticity moves about. Since under uniform rotation, m = 2$!, the time 
derivative of the curl of Eq. (9) is just this equation linearized. The longitudinal part of the left side of Eq. (9) is 
trivially the gradient of a scalar. Equation (9) constrains the number of degrees of freedom in two dimensions to four 
from the original five - n, w, and e. 

The time derivative of Eq. (9) is the superfluid acceleration equation, 

m + 20 X ij = -V(m - Kff), (12) 

where /x is the chemical potential. For an axially symmetric harmonic conflning trap of frequency wj, in the transverse 
direction and Wz in the axial direction, 

Kff=f [(c.i-^^')r^+u;,V], (13) 

where r denotes {x,y). In the frame corotating with the vortex lattice, the chemical potential fj, is related to the 
phase by 

MM)-K« = --4f^. (14) 
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B. Conservation laws 



The dynamics are specified by the conservation laws of particles and momentum, together with the superfluid 
acceleration equation, (12). The continuity equation takes its usual form, 

^ + V.J(M)=0, (15) 

where n is the density and j = nv the particle current. Conservation of momentum reads: 

a? - ^ ^ - 

+ 2mn X j + VP + nWeff = -&-(. (16) 

Here P is the pressure, and a is the elastic stress tensor, discussed below. At zero temperature, VP = nV/i, while in 

equilibrium, VP + nVVcs = 0. To calculate the displacement autocorrelation fimctions, wc include here, as in [23], 
an external driving force, — ^ (r, t), acting on the lattice, derived from an external perturbation H' = C, ■ e. 

The elastic stress tensor is derived from the elastic energy density of the lattice, which in two dimensions has the 
form (in the notation of [21]), 



f(r) =2Ci(V-e)2 + C2 







(dex_ 




V dx 


dy) ^ 


\dy 


dx J 



(17) 



where Ci is the compressional modulus, and C2 the shear modulus of the vortex lattice. The elastic constants are 
density-dependent properties of the fluid. Then the elastic stress tensor, cr,, is given in terms of the total elastic 
energy, Ee\ = J (fr£{r), by 

ai{r, t) = ^ = -4Vi (CiV • e ) - 2V • (CaVe,) , (18) 

where we allow for effects of density gradients entering through the elastic constants. In an incompressible fluid, 
C2 = nO/8 = — Ci. On the other hand, in the quantum Hall regime, the shear modulus is determined by the 
deviations of the interaction energy caused by distorting the vortex lattice, and [25] 

^ ^ms"n, (19) 

in agreement with the shear modulus given numerically in [11]. Calculation of the elastic constants Ci and C2 over 
the full range of f2 from the stiff Thomas- Fermi to the quantum Hall limits will be given in Ref. [25] . 

Subtraction of Eq. (16) divided by n, from the superfluid acceleration equation (12), with VP = nV/x, yields 

2ftx(k-v) = ^^. (20) 
mn 

Equations (12), (15), (18), and (20) fully specify the problem for a trapped system with a non-uniform density. 

C. Modes 

In the following, we neglect effects of non-uniformity of the equilibrium density for simplicity, and proceed to derive 
the modes as in [22] . The curl of Eq. (20) is 

while its divergence, together with (10), yields, 

^ x'e+2n^ ■e = -—^V-(a + C), (22) 
2Unm 
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where V x = -2C2V2(V x e), and V • ct = -2(C2 + 2Ci)V2(V • e). To eliminate v from (21), we note that the 
density oscihations are governed by 

-^+s2vMn = 2nOVxe, (23) 

where the sound speed, s, is given by ms^ = dP/dn [26]. In terms of the frequency, w, and wave vector, fc, we have 
then 

k-v= ^^'^ „,M -kx e, (24) 
so that in terms of longitudinal and transverse components [27] , 



C2 + 2Ci , 2\ , _ Ci 



iwer + 2f2 + -k-" Ul = - 



17r?,m J 2J7nm ' 

, 2a>2f^ C2 ,2\ Ct 



Solving for ez, and we have 



where the secular determinant, whose zeroes determine the mode frequencies, is 
D{k,uj)=(j^ -co^ 



4ii'+(s' + —{Ci+C2)]k' 
' nm ' 



+ ?f!^fc4 = (0,2 - ujj){uj^ - uj^) = 0; (27) 
nm 



we have dropped terms of second order in the elastic constants. 

For 2s'^C2k'^/nm -C (417^ + (s^ + 4(Ci + C2)/nm)k'^Y, as is always the case at long wavelengths in both the 
incompressible and quantum Hall limits, the mode frequencies are given by 



and 



^|=.4n2+(,2^ 4(Ci+C2) ,^2 (28) 
' nm ' 

2 ^ .2q^ 

'^'^ nm (4172 + (s2+4(Ci+C2)/nm)fc2)- ^ ' 

The first mode is the standard inertial mode of a rotating fluid; for Q, <C s^k"^ the mode is a sound wave, while 
for Vl ^ s^/c^, the mode frequencies begin essentially at twice the axial trapping frequency. This mode has been 
calculated in realistic trapping geometries in [28] and [29]. The second mode is the elliptically polarized Tkachenko 
mode. Equation (22) implies that the inertial mode is circularly polarized; thj^T — i, and in the Tkachenko mode, 
cl/ct — iivrf^^', the small longitudinally polarized component is 7r/2 out of phase with the transversely polarized 
component. In the limit of an incompressible fluid (s^ 00), 

uj = 4fl^+(s' + ^) k\ (30) 
\ nm J 

and the Tkachenko frequency, ujt, is linear in k, Eq. (7). In the soft limit, by contrast, 

2 s^C'2 4 , , 

unlike in the stiff Thomas-Fermi regime, the mode frequency is quadratic in k at long wavelengths; using Eq. (19) for 
C2 we have 

9 s^k^ 

ojt ^ ^ (32) 

The present results for the modes are valid for a uniform system over the entire range of rotation frequencies, from 
the slowly rotating stiff regime up to the melting of the vortex lattice. The Tkachenko mode has been calculated 
numerically for realistic trapping geometries in [19] in the stiff limit and more generally in [30]. 
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III. CORRELATION FUNCTIONS 



We turn to determining the effects of the lattice modes on the lattice ordering, and the phase coherence and 

condensate fraction of the rotating supcrfluid. To do so wc construct the correlation functions of the density, supcrfluid 
velocity, vortex displacements, and phase from the dynamical equations in the previous section. All the correlation 
functions of interest, including the single particle Green's function, can be written in terms of the density-density and 
displacement-displacement correlation functions. In the following wc let {AB){k, z) denote the Fourier transform in 
space and the analytic continuation of the Fourier transformation in imaginary time to complex frequency, z, of the 
correlation function {A{rt)B{r't')) - {A{H)){B{r't')). 

The density-density correlation fimction, {nn){k,z), is readily found from the response of {n{r,t)) to an external 
potential, U{r,t), coupled to the density. Using Eq. (24) to eliminate from (23), we have 

with D given by (27). As ^; — > oo, {nn){k,z) nk'^/mz'^, which is the expected f-sum rule on the spectral weight 

B{k,LO), 

/ ^B{k,^) = —. 34 
J-oo 27r m 

Similarly, {nn){k ^ 0,0) — * —n/ms^, yielding the correct compressibility sum rule, 

lim / ^-^^ = — t;. (35) 

fc^o 27r w ms^ 

As expected, the f-sum rule is dominated by the high frequency inertial mode; the low frequency Tkachenko mode 

dominates the compressibility sum rule [31]. 

The correlations of the longitudinal velocity are given in terms of {nn){k,u}), as usual, by 

2; 

{vLn){k,z) = —{nn){k,z), (36) 

and 

z2 1 

{vLVL){k, z) = ^{nn){k, z)-— (37) 
^ ^ z^(m-' + s\^ + —{2C^+C2)k^\-'^^k'\. (38) 



nmD \ \ nm j nm 

The correlation functions of the elastic displacements are given by 

du) Bij{k,Lu) _ S{ei} 



where (ci) is the displacement in the i^^ direction induced by the force (. Equations (39) and (25) then imply. 



nmD{k,z) \ 2Q?nm 

We note, for later calculation of the phase correlations, that the displacement-density correlation function, found 
from Eq. (23) together with the continuity equation, is 

(ner) {k,z) = {ern) {k,z)= ^J^^^^^^ {erer) , (43) 
and thus the displacement-longitudinal velocity correlation is, 

20^2 

{eTVL){k, z) = {vLeT){k, z) = ^2 _ ^2^2 {^tct). (44) 
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A. Lattice displacements 



We first address the effects of the vortex modes on the lattice displacements from equilibrium. At finite temperature, 
T, the equal time displacement correlations are given by 

((e,(r)-e,(r')f) = 2| -^0^ (^1 - cask ■ r) j^^ ^Bu{k,u;){l + 2f{co)), (45) 

where Z is the thickness of the system in the z direction, R = r — r and /(w) = l/(e^'^ — 1), with /? = l/ksT. The 
spectral weights Bij are found from Eqs. (40)-(42) by letting z ^ lj and 

The leading terms in the mean displacement of a single vortex from equilibrium due to excitations of the modes are 



d^k 1 



21.2 



+ 2/(0;,)) + 1^(1 + 2/(a;T)) 



(47) 



The mean displacement is convergent at zero temperature. However, at finite temperature it diverges logarithmically 
with system size if the Tkachenko mode spectrum reaches down into the soft quadratic regime; then 

The relative separation, with only the leading terms kept. 



((e(r)-e (/))') 



£k 1-cotik-R 
(27r)2Z nmujj 



2 1.2 



(49) 



converges at zero temperature, and lattice preserves long range positional order at large separation; at finite temper- 
ature, however, ((e(r) — e (r'))^) grows logarithmically with separation, ~ {Tm^/4:wZC2)inNy{R), where Nv{R) is 
the number of vortices within radius R. In the stiff Thomas-Fermi limit, this expression becomes (47r/nZA^) lnNy{R}, 
where A is the thermal wavelength; in the quantum Hall limit we have rather (207r^/81i^) (T/ms^) In A^„(i?). 

Equation (47) may be used with a Lindemann criterion to estimate the point where the lattice melts at T = in 
an extended system [10,11]. The zero temperature displacements are sensitive to the entire spectrum of modes up to 
the lattice Debye vector, kj, = (47rn„)^/^, whereas the mode frequencies derived here are valid only for k <C kd- For 
a first estimate, we replace the integrand by its infrared limit, letting / (Pk/{2w)^ riy. Then using the Tkachenko 
frequency at rapid rotation, (32), we have 

/?2\ 1 / „2i,2 

(50) 



il2 = L(i~^\ 

P 2v\ AujT^ ) 




where (. = (I/tto^)^/^ is the radius of the Wigner-Seitz cell around a given vortex, and we have used (19). The 
"1" term arises from the inertial mode, with equal contributions from the transverse and longitudinal displacements; 
the final term arises solely from the soft Tkachenko mode contribution to the transverse displacements. In order to 
take into account approximately the mode structure at larger wave vector, we include the kinetic term — V^n/4m 
in the pressure. In the non-rotating weakly interacting gas, this term modifies the linear spectrum, sk, into the 

1 /2 

full Bogoliubov spectrum, Ek = ((sfc)^ + (A:^/2m)^) . Inclusion of such a term is equivalent to replacing by 
+ k'^ /Arn^ in the mode spectrum, with the effect of stiffening the Tkachenko mode at high wave vector. Then 
at fl/ms"^ ^ 40, as expected under typical experimental conditions at ;/ ^ 10, the contribution of the Tkachenko 
modes to (e^) is reduced from 1.73/i^ in Eq. (50) to ^ 0.72/i^, a result consistent with that reported in Ref. [11], 
mfl{e^} ~ 0.66/j/. Then with the contribution of the inertial modes included, 

(e^) 1.22 

Taking the Lindemann criterion for melting in two dimensions put forth in [10], (e^)/^^ ~ 0.07, we find melting of 
the lattice at filling i/ ~ 17; were we to take only the Tkachenko mode contribution, the melting would be at ^ 10. 
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B. Phcise correlations 



We next determine the effect of the vortex excitations on the correlations of the order parameter in the superfluid 

and on the condensate fraction. The condensate density, no, is most conveniently found, a la Onsager-Penrose, as the 
limit of the single particle density matrix for large separation: 



(V(r)V^(r')) 



no, 



r — r 



oo, 



(52) 



where ipif) is the single particle annihilation operator [32]. We write ip(r) in terms of the density and phase operators, 
viz., ^(r) = y^n(r)e**(''\ and expand the long wavelength structure to second order in small fluctuations of n about 
n, its equilibrium value, and e**^'"^ about unity. Then in terms of equal time correlation functions. 



(V'(r)V^(r')) ^ n(e**('^)e-**('^')) + ^ {{5n{r)5n{r')) - {{5n{r)f)) 



(53) 



where 5n{r) = n{r) — n, and n denotes the average density. The first term on the right is the U(l)-invariant correlation 
of the order parameter, given, for Gaussianly-distributed Fourier components of the phase, by 



(e 



(54) 



Equation (33) impHes that as |r — r '| — > oo, {6n{r)5n{r')) and also the final bracket in (53) vanish. Thus the density 
of particles in the condensate is given by. 



no = lim 

|r — r'l — ^oo 



ne 2 



i<(3.(f)-$(f')) = ) _ ± 



4n 



{Snirf) 



(55) 



When the phase flucutations are convergent, and vanishing at large separation, we can expand to second order to find 
the usual expression for the condensate depletion in the Bogoliubov approximation (see Sec. IV): 



n' = n{^'} + -{Snirf}. 



(56) 



To determine the phase-phase and the density-phase correlations, we note that the divergence of Eq. (9) implies 
that the Fourier compenents of the phase obey $ = —{im/k){vL — 2fleT)- Thus 



m 



[{vLVL){k, z) - 2il {{eTVL){k, z) + {vLeT){k, z)) + An''{eTeT){k, z)] 



which, with Eqs. (38), (41), and (44), becomes. 



($$)(fc,^) = ^ 



n^fc^ 



[nn) 



nk 



4f22 



+ s^k"^ 
— s^k'^ 



nD 



■402 



4(Cl+C2)^2 



nm 



Similarly, 



{^n){k,z) = -i^ 



{nn){k, z) 



z^ — s^k"^ 



-iz — ^{<^^){k,z). 



(57) 



(58) 



(59) 



Note that the phase fluctuations are more divergent in the infrared limit by a factor l/fc^ than the transverse dis- 
placement fluctuations. 

The relative phase correlations are given by: 



d^fc 1 — cosfc-.R 



(1 + 2/(u;t)), 



(60) 
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plus finite terms, where R = r — . In the stiff limit (s^fc^ ^ il^), the relative phase fluctuations are convergent at 
zero temperature, but, as discussed in [23], the finite temperature contribution of the Tkachenko modes washes out 
the phase correlations logarithmically, and the system is in fact in a fragmented condensate phase [33]. One finds 
similar behavior at zero temperature in the soft limit. The zero temperature part of the integral (60) is ^ ln{kDR)- 
Thus from Eq. (54), we find the correlation of the phase factors, 

(e**We-**("')) ~ {kDR)-" ~ KiR)-"'^, (61) 
where again Ny{R) is the number of vortices within radius R, and 

/ o_x 1/2 

Prom the limits Q. <C ms^ toVt^ ms^, the range of ry is, 

>v> — ; (63) 



u \ n J - ' - 

is . The phase correlations fall off algebraically at large i?, as expected for a two-dimensional system [34]. As the 
phase correlations decrease, the condensate fraction also falls, as (rio/n) ^ Ny"^^"^ . 

The falloff of the phase correlations begins to become important for {rj/2)\nNy > 1, which translates effectively 
into the condition in the quantum Hall limit that v ^ IJlnNy, or N ^ z/e"'^''; for A'' = 10^, the condition is that 
Nyk,5x 10^; for N = lO", > 10^; and for N = 10^, one needs only Ny > 10^ to find loss of phase coherence 
across the system. 

One may ask whether the falloff is significant by the time the lattice will have melted. The most divergent terms 
in the phase fluctuations are induced by the transverse displacement fluctuations in Eq. (58), so that in a system of 
transverse radius R±, one has roughly, 

\rJ^n^ - ^ |(^"'>(4), (64) 

where the average {l/Awny) InNy is taken with respect to the weight in the transverse displacement correlation 

function. Thus 

lim i(($(r)-$(r'))') ~2^1n7V,. (65) 



r — r' 



The Lindemann criterion, (e|,)/€^ ~ 0.07, implies that for Ny ^ 10^ the right side of Eq. (65) exceeds unity at the 
melting point, setting the scale for number of vortices for which loss of long range order of the condensate prior to 
melting becomes important. 

At finite temperature, the phase correlation integral (60) is logarithmically singular in the infrared for all R, 
indicating that in two dimensions, the system is, as expected, no longer Bose condensed. The situation is the same as 
in the stiff Thomas- Fermi limit, where the phase correlations in two dimensions diverge in the infrared, but in three 
dimensions fall algebraically with separation [23] . 



IV. THE SINGLE PARTICLE GREEN'S FUNCTION AND CONDENSATE DEPLETION 

We now determine the structure of the single particle excitations in terms of the modes of the lattice by constructing 
the long wavelength behavior of the single particle Green's function, 

G{rt,r't') = -i{T[{iP{rt) - {ip{rt))){tp^r't') - (V'^(r-Y))]) (66) 

(T denotes the time ordered product), in terms of the correlation functions calculated in the previous section. Again 
expanding in small fluctuations about equilibrium, we have 

G(k, z) = ^(nn)(k, z) + ^ ({ne-'^)(k, z) + (e'*n)(fc, z)) + n(e'*e-'*)(fc, z). (67) 
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To second order in the fluctuations of the density and phase, 

G{k,uj) = ^{nn){k,z) + ^{{^n)ik,z)-{n^)ik,z)) + n{^^){k,z). (68) 

To illustrate this method of calculating G, we first consider a weakly interacting non-rotating system, for which 
Eqs. (33), (27), (59), and (58) imply. 



and 



{nn){k,z) = -^^, (69) 

777 2; 

($n)(fc, z) = -(n$)(fc, z) = -i—{nn){k, z), (70) 

7D ^ / 771 7 \ ^ 771 

($$)(fc,z) = —iv,v,) = (-^) (nn)(fc,z) - ^, (71) 

1 /2 

where Ek = (^gnk'^/in+ (fc^/2TO)^) is the Bogoliubov single particle energy. (The term ^ k^ in E"^ appears only 
if one includes the kinetic term, —V^n/4m in the pressure.) Substitution of these correlation functions into Eq. (68) 
yields the usual Bogoliubov result: 

G{k, z) = z-^-El ■ ^ ^ 

In the presence of rotation, we use Eqs. (58) and (59), to write 

Gik, = + „..^) - 4n^ - 4£i±^) + |1 - 2^^^) I . ,73) 

D [ V nm)Am\ nm J } 

The contribution of the modes to the density, n^, of particles of momentum k excited out of the condensate, to leading 
orders in fc^, is then 



LUT \ 2 / SmO, 



,'ms'^\\ ... If ms^ 



(74) 



The contribution of the inertial mode (the final terms) is always finite in two dimensions. However, the contribution 
of the Tkachenko mode in the soft limit (the first term) leads to a logarithmic divergence of n' in two dimensions 
at zero temperature [11], and a quadratic divergence at finite temperature, thus destroying the condensation in an 
infinite system. 

It is instructive, finally, to calculate the effect of the lattice modes on the superfiuid mass density, ps, which is 
related to the transverse velocity autocorrelation function by [35] 

Pn = P - Ps = -{mnf \mi{vTVT){k,LO = 0). (75) 
fc— >o 

Since vt = —2Q,eL, we find from Eq. (40) that 

Pn = mn, (76) 

that is, the superfiuid mass density vanishes. As discussed in [23] the lattice excitations replenish the sum rule (75), 
a reflection of the fact that the moment of inertia of a rapidly rotating superfiuid is effectively the classical value [36] . 
The vanishing of ps is consistent with the behavior of G{k 0,0), which according to Josephson's sum rule should 
approach —m^rio/ p^k"^ [35], while in fact G{k — > 0,0) — > —{Un/C2){2w?/k^). 
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